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The possibility of a significant effect of exotic particles on the Z' — > 777 and Z — > 777 decays is 
investigated in the context of the minimal 331 model. This model, which is based in the S(7c(3) x 
SUl{3) X Ux{l) gauge group, predicts the existence of many exotic charged particles that can 
significantly enhance the decay widths. It is found that the standard model prediction for the 

■ Z — >■ 777 decay remains essentially unchanged, as the new physics effects quickly decouples. On the 
other hand, it is found that the contributions of the new exotic quarks and gauge bosons predicted 

■ by this model lead to a branching fraction for the Z' — > 777 decay of about 10^^, which is about 3 
^Nj ■ orders of magnitude larger than that of the Z — >■ 777 decay. 
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I. INTRODUCTION 



Numerous extensions of the Standard Model (SM) predict the existence of new neutral gauge bosons Z' 1]. This 
_ , class of gauge bosons can be associated with spontaneous symmetry breaking (SSB) of additional U{l)s gauge groups or 
with Kaluza-Klein excitations of theories with extra compact dimensions 0] ■ Phenomenologically, the most interesting 
' option is the breaking of these U{l)s at around TeV scales, giving rise to extra neutral gauge bosons observable at 
1^ . the Large Hadron Collider (LHC). Although the LHC has been designed to solve the Higgs puzzle, it is possible that 
in the first stage of its running signals of new physics show up, such as the resonance produced by a new Z' gauge 
boson decaying into charged leptons Therefore, it is worth investigating the phenomenology of this particle. In 
CN ' this paper, we will focus on the Z' gauge boson predicted by the minimal 331 model which predicts new physics 

^ , at the TeV scale. In particular, we are interested in studying the rare decay of this Z' gauge boson into three photons. 
• These types of decays are naturally suppressed in renormalizable theories, as they first arise at the one- loop level. 
The analogous decay in the SM Z 777 is quite suppressed, with a decay width of the order of 10^^° GeV p-[a]- 
the SM, the fermionic contribution was calculated in an approximate way more than two decades ago \d\ and some 
years after the complete analytical results were presented [7|. As the W gauge boson contribution is concerned, it was 
presented in 0, It was shown in Q that the fermionic- ly interference is significant. Since there are no charged 
scalars particles in the SM this kind of contribution only occurs in models with extended Higgs sectors. In P such 
a contribution is analyzed in the context of the popular two Higgs doublet model and the minimal supersymmetric 
standard model. It was found that this contribution is several orders of magnitude lower than the fermionic and W 
contributions. A more general analysis of the contribution of all three types of particles to this decay was presented 

5^ ■ As already commented, we are interested in studying the Z' — > 777 decay in the context of the minimal 331 
. model [1,1^. This model, which is based in the SUc{i) x SUl{3) x Uxi'i^) gauge group, has attracted the attention 
of numerous authors in the past decade, mainly because it possesses some peculiar features that are not present in 
other SM extensions. Its more interesting property is that anomalies do not cancel in each generation independently 
as in the SM, but such a cancelation occurs only when the three generations are taken into account together. The 
cancellation mechanism requires that the number of families must be an integer multiple of the color number. This 
fact, together with the property of asymptotic freedom of QCD, which establishes that the color number is less than 
five, implies that the 331 model predicts the existence of only three fermionic families. Another interesting property 
of the minimal version of this model is that the weak angle is subject to the constraint < j It results that, 
when it is evolved to high values, the model loses its perturbative character at a scale about 8 TeV The fact that 
the value of is very close to 1/4 leads to an upper bound on the scale associated with the first stage of SSB, when 
the SUci^) X SUl{3) x Ux{l) group is broken down into the SM group SUc{3) x SUl{2) x C/y(l), which translates 
directly into the gauge bosons that acquire masses at this scale, among them the Z' gauge boson tSjiilS]- Then, the 



^ From now on, and cyy stand for sine and cosine of the weak angle 6y 
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331 model is phenomenologically well motivated to be probed in the LHC. 

In the minimal 331 model the lepton spectrum is the same as in the SM, but it is arranged in antitriplets of 
the gauge group SUl{5)- The quark sector is also arranged in the fundamental representation of this group, which 
requires the introduction of three new quarks. Along with the Z' boson, four additional charged gauge bosons are 
predicted by the minimal 331 model: two singly charged bosons and two doubly charged ones Y^^. These 
gauge bosons carry two units of lepton number and so have been classified as bileptons [l^ ■ The new gauge bosons 
together with the exotic quarks are endowed with mass at the first stage of SSB, when SUl{3) x [/x(1) is broken 
into SUl{2) X t/y(l) [l|. Since SUl{2) is completely embedded into 5'[/l(3), the couplings between the SM and 
the extra gauge b osons are determined by the coupling constant g associated with the SUl{'2.) group and the weak 
angle Ow [iJ, ^M- On the other hand, the Higgs sector of the minimal 331 model consists of three triplet and one 
sextet, but only one triplet is needed to break down SUl{2>) x Ux{^) into SUl{2) x Uy{^)- At this stage, the Higgs 
sector consists of three doublets and one triplet of SUl{2), as well as of four complex singlets [l^,!!^- After breaking 
the usual electroweak group into the electromagnetic one, the physical scalar sector is composed by: five neutral 
CP-even Hi{i — 1, • • • ,5), three neutral CP-odd Ai {i — 1, 2, 3), four charged hf (i ~ 1, • • • ,4), and three doubly 
charged df^ {i — 1, 2, 3), from which only Hi is light, with mass of the order of the Fermi scale v [T^. In a previous 
communication by some of us [l7j, a comprehensive analysis of tree- level two-body decays of Z', including the one- loop 
induced ones Z' — > Zj and Z' — > ZZ, was performed in the context of the minimal 331 model. Our main goal in 
this work is to investigate the sensitivity of the F° 777 {V^ — Z, Z') decays to virtual effects of new particles of 
spin 0, 1/2, and 1 living at the TeV scale. As already commented, one interesting peculiarity of the 331 model is the 
presence of new particles with a charge content that differs from that carried by the known particles. Exotic particles 
with a charge content Q > 1 in units of the positron charge may enhance substantially the amplitude of the ^^777 
vertex, as it is proportional to Q^, which in turns leads to a probability proportional to Q^. This class of effect was 
studied by some of us some years ago in light by light scattering, showing that this process is quite sensitive to the 
contribution of the doubly charged bileptons, as they modify the cross section by an additional factor of = 2^ p^ . 
In some more recent communications by some of us, the impact of the exotic quarks on the rare ^ 999 [13 ^'^'i 
~^ 991 decays was studied. All the new charged particles that predict the model receive their mass at the first 
stage of SSB, so that a good strategy to estimate the impact of scalars, fermions, and vectors to the — ?> 777 decay 
is to assume that each set of these three types of particles are mass degenerate. Accordingly, the impact of the three 
exotic quarks, two of them {D and S) with charge content of Qd = Qs = —4/3, and the third one (T) with charge 
Qt = 5/3, is + Q| -I- 32.7. On the other hand, the impact of the simple and doubly charged bileptons 

is Qy±± + Qy± = 65. As far as the four simply charged and three doubly charged Higgs bosons is concerned, their 
impact would be proportional to 3Q^±± + 4(5^± = 196. 

The rest of the paper has been organized as follows. In Sec. |ll]a brief discussion of the 331 model is presented. 
Section IIIII is devoted to present compact analytical expressions for the amplitudes associated with the — > 777 
decay, with V'^ = Z, Z'. In Sec. IIV[ we discuss our results. Finally, in Sec. |V]the conclusions are presented. 



II. THE MINIMAL 331 MODEL 



As already mentioned in the introduction, the lepton spectrum of the model is the same as in the SM and are 
accommodated as antitriplets of 5'C/z,(3): 
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Notice that (i?^) ^ = {iaR)'^ (^o — e, ^, r). In order to cancel the S'[/i,(3) anomaly, the same number of fermion triplets 
and antitriplets is necessary. This requires one to arrange two quark generations as triplets and the other one as an 
antitriplet. It is customary to choose the third generation as the one transforming as antitriplet in order to distinguish 
the new dynamic effects in the physics of the quark top from that of the lighter generations. Accordingly, the three 
generations are specified as follows: 
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Q3^ t , (3, 3*, 2/3). 

On the other hand, the right-handed quarks are specified as fohows: 

d^ : (3*, 1,1/3), : (3*, 1,4/3), 



(3) 



(4) 



^i^ c", t" : (3*, 1,-2/3), : (3*, 1,-5/3), 



(5) 



In the first stage of SSB, when the S'C/l(3) x Uxi^) group is broken into the usual electroweak group SUl{2) x C/f(1), 
only the three new quarks D, S, and T acquire masses. These exotic quarks arise as singlets of the SUl( 2) g roup, so 
they do not couple to the W gauge boson. However, they do couple to both the Z and Z' gauge bosons [l7| . 

On the other hand, as already commented in the introduction, the Higgs sector of the minimal 331 model is 
comprised of three triplets and one sextet of SUl{3): 
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To break SUl{3) x Ux{^) into SUl{2) x J7y(l), only the scalar triplet of SUl{3) is required. The hypercharge is 
identified as a linear combination of the broken generators and X: Y = ■\/3(A^ + V^XX^), with a Gell-Mann 
matrix and = -y2/3diag(l, 1,1). The next stage of SSB occurs at the Fermi scale and is achieved by the two 
triplets and 02. The sextet H is necessary to provide realistic masses for the leptons [llj. In these expressions 
, 'i'l , ^2 = icr^$2 ^^'^ ^3 s-^^ all doublets of SUl{'2.) with hypercharge 3, 1, 1, and 1, respectively. On the other 
hand, T is an SUl{'2) triplet with Y = +2, whereas S~ , p , and -q are all singlets of SUl{^) with hypercharge —2 
—4, and +4, respectively [iJl- The extra Z' boson, the bileptons and the exotic quarks get masses at the first stage 
of SSB through the vacuum expectation value < (j)Y >o= (0, 0,m/\/2). The bileptons form a SUl{2) doublet with 
hypercharge +3. The spectrum of physical gauge particles is the following. The charged gauge bosons are given by 



Y++ 



Y+ 



1 



^(4-^4) 



(8) 
(9) 



W+ 



V2 



(4-^4): 



(10) 



with mY++ 



g2(u2+i'| + 3u|)/4, 



Y+ 



g'^{u'^ + vl+vl)/4, and mp^, = g'^ {vf + V2 + v'^) / A. The hierarchy of the SSB 

< 3m^. However, to simplify the discussion 
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yields a splitting between the bilepton masses given by |my+ 
we will consider only the degenerate case. 

In the neutral sector, the gauge fields {A'^, A^, X) define three mass eigenstates {A, Zi, Z2) via the following rotation: 
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where the mixing angle is 

9 2 

sin^ 9 = |i- , (13) 

m^^ — m^^ 

with = ''TT'w/'^w ^^^'^ -^1 standing for the SM Z boson. In this paper, we will work in the approximation 9 — 0, 
so Zi and Z2 coincide with Z and Z', respectively. 

To calculate the amplitudes for the — >■ 777 decay, we need the Feynman rules for the couplings of Z and Z' to 
all the charged particles of the model. The couplings of Z and Z' to leptons and quarks, including the exotic ones, are 
all given in 17] and we refrain from including them here. The couplings of these particles with photons are dictated 
by spinorial electrodynamics. As far as the couplings of the neutral gauge bosons A, Z, and Z' with the charged ones 
W^, F+, and Y~^^ , they depend on the gauge-fixing procedure used to quantize the theory. The calculation of these 
contributions are greatly simplified if one uses a nonlinear gauge- fixing procedure. To carry out the contribution 
we used a covariant gauge- fixing procedure as the one presented in [isl [22j . As far as the bilepton contribution is 
concerned, we used the Feynman rules that arise from the gauge- fixing procedure used in (23l, [2^. The Feynman rules 
for all possible couplings among the neutral and charged gauge bosons of the minimal 331 model can be found at the 
Appendix [A] Some of these gauge bosons rules have been worked in (isl [TtI . [23| . On the other hand, the coupling 
of charged scalars to photons are model independent, as they are dictated by scalar electrodynamics. Consequently, 
we only need the couplings of all the charged scalars with the Z and Z' gauge bosons. These couplings arise from 
the Higgs kinetic sectors of the 331 model, but to determine them one needs to diagonalize the Higgs potential of 
the model. We have used the diagonalization given in [l^ to determine all couplings of the Z and Z' gauge bosons 
with the charged scalars of the model. The vertex functions associated with the V^°"S't(fci)5(fc2), V^"A^S^S, and 
A^A^S^S couplings, with S stands for a charged scalar, are given, respectively, by 

^7-^(K1 - k2)a , gap , «2e Qsffa/? , (14) 

ZCw cw 

where all momentum are taken incoming to the vertex. The values of and gy are presented in Table ID In this 
table, 9 is the Z — Z' mixing angle, which we have taken equal to zero. On the other hand, the H22, H32, a, N4, N^, 
X4, and are given by 

^^22- , H32^--~ , a=V2/ , (15) 



V2 ' f 
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32/2-ka(at\/32/2 + A2 



^4.5 = -J(at\/32/2 + A2) , (16) 



where A and / represent parameters of the Higgs potential 16], which we have assumed of the same order: A = Ag = 
Ai7 = Ai9 - 0(1) and / = /i = - 0(1). 

III. THE AMPLITUDE FOR THE V° 777 DECAY 

In this section, we present the amplitudes for the on-shell vertices ^^777, with V'^ — Z,Z'. Our notation and 
conventions are established in Fig. [1] These couplings first arise at the one-loop level through diagrams shown in Fig. 
[2I in which circulate all charged particles of the model. The invariant amplitude can be written as follows: 

■^777V» = E-^7™'''^A.ibi- Ai)e;.,(p2,A2)e^3b3^ , (17) 

X 
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TABLE I: Values of Qs and gyo in the couplings of Z and Z' to charged scalars. 
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where X denotes the type of particles circulating in the loops. The fermions will be collectively denoted by F = u, d, 
s, c, 6, t, D, S, T, e, fi, t. On the other hand, gauge bosons and their associated pseudo Goldstone bosons and ghosts 
contributions will be separated into vector and scalar contributions, namely, vector {V) and scalar [S) particles, with 
V = W+,Y+, y++, and 5 = 6+ , C+ , (7+ , G+, C+, (5+, G++, C++, C++. The contribution of the physical Higgs 
bosons will be denoted hy H — h'^ , , , , , 4^- ^'^'^ calculation purposes, it is convenient to organize 

the amplitudes according to the spin of particles circulating in the loops. This is possible since we calculate the gauge 
particles contribution using covariant gauges -R^-gauges [l5l [22l - [2^ . which separately render finite and gauge invariant 
the pseudo Goldstone bosons and ghosts contributions 

M^^^v" =^Mi+Mi+Mo , (18) 

where A4i, Mi, and A4o are, respectively, the spinorial, vectorial, and scalar amplitudes, which can be written as 
follows: 

Mi^^Mf, (19) 

F 
V 

Mo^Y.-^s,H- (21) 

Structurally speaking, these are the only amplitudes that are different. The Lorentz tensor structure of the amplitudes, 
as well as the functional way of the form factors involved, are dictated by electromagnetic gauge invariance and Bose 
symmetry. Gauge invariance requires that 

P.M. E K^^v''"^" = , z = 1, 2, 3 , (22) 

X 

whereas Bose statistics dictates that the amplitudes must be symmetric under the interchanges of pairs of photons: 

(Pl,A*l) (P2,Ai2) ^ (P3,M3) • (23) 

There are six different configurations by each type of diagram. 

The fermionic contribution is given by six box diagrams as that shown in Fig. [Ifa). The corresponding tensorial 
amplitude is given by 

^MiM2M3M4 ^ J2 Wp^'^bZ'^' . (24) 
^ F 
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In the case of the vector contribution all types of diagrams in Fig. [2] contribute. The tensorial amplitude can be 
written as 

A^MlA'2M3A»4 _ I A^MlP2/J3A'4 , A^Ml/J2/i3M4 , . ^ /il A12A13 AI4 , A^A»lAi2M3A'4 ) (r,r\ 

■'^'1 — / . y^'-V.Box ■^^'■V,Trianglel + ■^^W,Triangle2 ^ ■^^'■V,Bubble J ■ 1^°/ 



The Feynman rules needed for the vector contribution are given in Refs. |15l . I22l - l24| , these and other new rules are 
summarized at the Appendix [A] 

Finally, the scalar contributions are characterized by the following amplitudes 

A/<MlM2M3M4 _ ( .^A'1A'2M3M4 I .^AtlAt2A»3A»4 I A/<A'1M2M3A'4 I .^MlM2M3M4 ) /Of?\ 

'0 — 2^ [■^^^X,Box ■^^^X,Trianglel ^ ■^^^X,Tnangle2 ■^^'■X, Bubble j ' 

X ^ ^ 

here X stands for a nonphysical scalar S' or a Higgs boson H . The diverse vertices involved in the above amplitudes 
were given in the previous section for the case of the couplings of with scalars, whereas the corresponding couplings 
of the photon are dictated by scalar electrodynamics. 

Once the loop integrals are solved, the amplitudes can be expressed in terms of gauge structures and their associated 
form factors as follows: 
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i—l 

where 

FyOj = F-j^o^ + FyOi + FyOi . (28) 

The Lorentz tensors x^^^^^^^^ are gauge structures given by [l^ 

yP,P2M3M4 ^ . p^^,.U<2 _ p^lpf . p3gA<3P4 „ pP3pM4) ^ (29) 

J,|.U.2M3M4 ^ . p3^p,p2 _ pMlpA^2)(p^ . p3^WP4 _ pM3p^4) 

+ • P2Pr - Pi • P3Pr - Pfs^^^^) • (31) 

Using Bose symmetry, each of these tensors determines a set of 6 gauge structures and their associated form factors: 

r p rTiA»lA'2M3M4 p ^^MlA'2M3^'4 \ 

r T? TiMlA»2A'3M4 p TiMlA«2A'3A'4 \ 

■[J'yOyJy , i^y0i2J 12 /, 

As it can easily be verified, all these Lorentz tensors are gauge structures, as they satisfy the transversality condition: 

yMiM2M3P4 ^ , i = l,2,3. (33) 

Notice that the amplitude (P?]) is also transverse with respect to the gauge boson, which means that it appears in 
the T^°777 coupHng not directly, but only through the strength tensor V^^, — d^V^ — d^V^. 
The fermionic contribution is given by 
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A/<PlP2A'3P4 _ * 9^ \ ^ p2 Tn/il/i2M3M4 (OA\ 

■^l -^^l^^vo.T. , (34) 



=1 
where 

F=q,Q,l 

= Fyot + ^yoj + Fyo^ , (35) 

with the definition 

gl = -NpQlgU , (36) 
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being Np the color number, 1 for leptons and 3 for quarks, Qf is the charge content in units of e, and gyo is the 
vector-hke couphng of to fermions, which is given in In ([55]) q, Q, and I stand for known quarks, exotic 

quarks, and charge leptons, respectively. The form factors /yoj are exactly the ones given in [l9l | and we refrain from 
present here. 

As far as the gauge boson contributions is concerned, we first note that in nonlinear gauges the contributions of 
pseudo Goldstone bosons, ghosts and antighosts satisfy the following relations 

Mg = -Mc = -Mc ■ (37) 
Then, the spin 1 contribution comprises the vectorial amplitudes plus the nonphysical scalar ones: 

7Wi=;^Xv,s. (38) 
v,s 

Once the corresponding vectorial and scalar amplitudes are collected, the tensorial amplitude can be written as: 

3 18 

■^1 -^^Z.^vo,T, , (39) 
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where 



^V°t — X] 9l fv°i ^ X] 9ofvoi 

V=W+,Y+,Y++ S=Gl-.,G+,G+ + 

= F^o, + F^o,, (40) 



with 



gY = Qlg^o , gl = Q%g^o . (41) 



More explicitly, 



^Vi — '-^w+ \9vo Jv'H ~ 9v« Jvoi ) ' I'^^y 



Fv°i ~ {Qy+9v° + Qy++9vo — 

{qU9vI' + QU+9v^^) f% ' (43) 

and niY = mY+ ~ mY++ ■ Particularly, 

= '2cl,fl-C2wf'ir , (44) 

Ff,. = (45) 

Fj, = (7 - 34^2^) (/J, -/f-) , (46) 



Fl, - 9%/3 Jl - ^sl^fl., - sVScwf^,^, ■ (47) 



The various form factors appearing in the above expressions are given in the Appendix [ 
On the other hand, the contribution of the charged Higgs bosons is given by 

^yoi — 2^ 9q Jv-h 

H=hf Ji^ ,hi ,df + ,di + ,dj + 



F^o, , (48) 
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with 



where 



9o = Qh9v° 



(49) 
(50) 



H _ fh 

and rriH = iTif^+ — m^++. The form factors /^^ and are given in the Appendix IB] 
The differential decay width is given by 



(51) 



\M\' 



dxdy 2567r3 

where x = 2p\lmyo^ y = 2p2/raya^ z = 2p^/rayo, with x + y + z = 2. This parametrization leads to pi ■ p2 
rriyo {x + y — l)/2 , pi ■ ps — rriyo (1 — y)/2 , P2 ■ Pz — "lyo (1 — x)/2 , so the allowed region is determined by 1 
limits < a; < 1 and 1 — x < y < 1. This leads to a decay width given by 



(52) 



(53) 



where 1/3! is a symmetry factor due to the presence of three identical particles in the final state. In the above 
expression, 
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Vx^Y. FvhT^h''''^"^1. (Pi' Ai)e;,(p2, A2)e;3 (pi, A3)e^,(pi, A4 



(54) 



(55) 



with X denoting the virtual particles contributing to the — > 777 transition. The squared of the amplitude takes 
the form 
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xe^a (P3, A3)e^3 (P3, A3)e^^ {p4, A4)e*^ (p4, A4) , 



(57) 



where the factor 1/3 results from averaging on the polarization states of V^. In terms of the Vx amplitudes, the 
decay width can be written as follows; 
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FIG. 1: Notation and conventions for the V''777 vertex. 




FIG. 2: Generic Feynman diagrams contributing to the — >■ 777 decay. Fermions only contribute through box diagrams as 
shown in (a). 



IV. RESULTS AND DISCUSSION 



We now turn to discuss our results. We first discuss the impact of new physics on the Z 777 decay. We have 
verified that any virtual contribution X such that mx > mz/'i is marginal. To estimate these new physics effects, we 
consider a scenario with mg^SOO GeV and mY='mH='250 GeV. The relative importance of each type of contribution 
is shown in Tables and IIIII Our results for the SM contributions is in perfect agreement with those given in the 
literature [6-8]. From these tables, it can be appreciated that the contributions of the u and c quarks dominate. The 
contributions of the W gauge boson, as well as those arising from the charged leptons and the d, s and b quarks are 
all of the same order of magnitude and one order of magnitude lower than those given by the u and c quarks. It can 
be seen that the top quark contribution is marginal indeed, as it is two orders of magnitude lower than those induced 
by the u and c quarks. As far as new physics effects are concerned, the effect on this decay is quite suppressed, as 
their contribution is three orders of magnitude lower than the SM contribution, at the best. 

We now turn to discuss our results for the Z' — >■ 777 decay. Before analyzing the diverse type of contributions to 
the branching fraction of this decay, it is interesting to examine the role played by the spin of the particles circulating 
in the loop. As it was discussed in the previous section, the diverse contributions to the Z' — 777 decay can be 
grouped in accordance with the spin of the particles circulating in the loop, namely, spin (charged scalars), spin 1/2 
(charged leptons and quarks), and spin 1 (simple and doubly charged bileptons). Then, we analyze the behavior of 
the squared amplitude for only one type of particle as a function of the mass ratio m/mz', with m the mass of the 
particle circulating in the loop. The squared amplitude which we will examine is given by 



|M|^ 



(59) 
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TABLE II: Contributions to the Br(Z — >■ 777) in the scenario mg = 500 GeV, my — mn = 250 GeV. Interference effects 
among sectors are also shown. 



Sector 


Br 


Fermions 


4.16 X 10"^° 


Gauge Bosons 


1.03 X 10~" 


Scalar 


3.04 X 10"^^ 


Fermions- Gauge Bosons 


9.92 X 10"" 


Fermions-Scalar 


-5.90 X 10"" 


Gauge Bosons-Scalar 


4.30 X 10"^^ 


Total 


5.26 X 10"^° 



TABLE III: Fermionic and bosonic contributions to the Br(Z — )■ 777) in the minimal 331 model. Some interference effects are 
also shown. 



Fermions 


Br 


Quarks 


Br 


SM Quarks 


Br 


Bosons 


Br 


Quarks 


2.67 X 10~^° 


SM quarks 


2.67 X 10"^" 


u 


3.76 X 10"" 


W boson 


1.03 X 10"" 


Leptons 


1.66 X 10~" 


Exotic quarks 


6.70 X 10~" 


c 


4.08 X 10"" 


Bileptons 


1.88 X 10"" 


Quarks- Leptons 


1.33 X 10""' 


SM-Exotic quarks 


-3.62 X 10"" 


t 


1.11 X 10""'* 


W boson-Bileptons 


-1.57 X 10"" 










d 


1.91 X 10"^^ 














s 


1.91 X 10"^^ 














b 


2.12 X 10"^^ 














Interference 


1.82 X 10""' 






Total 


4.16 X 10"^" 




2.67 X 10"^" 




2.67 X 10""' 




1.03 X 10"" 



where 
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In this expression, /^Fj'" are numerical factors, which are irrelevant for the present discussion. In Fig. [3J the behavior 
of |Mp as a function of m/mz' for the three types of spins is shown. From these figures, it can be appreciated that in 
the case of spin 1/2 the highest contribution occurs for mf /mz' = 0.03, which is about of 1.4. This means that the 
mean contribution from the fermionic sector to the Z' — >■ 777 decay arises from the lightest SM charged leptons and 
quarks. On the other hand, from the second graphic of the same figure, it can be seen that the highest spin effect 
(0.15) occurs for niH/mz' — 0.12. This shows that the main contribution would arise from a relatively light charged 
Higgs boson, with a mass of about a 12% of the mz' mass. It is interesting to see that the scalar contribution is about 
of one order of magnitude lower than the fermionic one. As far as the spin 1 contribution is concerned, it deserves 
special attention. In first place is the fact that the Z' particle does not couple directly to pairs of SM W gauge bosons 
at the tree level, but only very weakly through the Z' — Z mixing. The absence of a direct coupling Z'WW is a 
consequence of the fact that the Z' gauge boson emerges in the first stage of spontaneous symmetry breaking as a 
singlet of SU{2)l. So, the only spin 1 contributions to the Z' — 777 decay arise from the simple and doubly charged 
bileptons. On this matter, one interesting feature of the minimal 331 model is that the new gauge boson masses are 
bounded from above due to the theoretical constraint which yields s^,/ ^ 1/4 [l2| • The fact that the value of is 
very close to 1/4 at the mz' scale leads to an upper bound on the scale associated with the first stage of SSB, which 
translates directly into a bound on the Z' mass given by mz' ^3.1 TeV [12], which in turns implies that the bilepton 
masses cannot be heavier than mz' /2 
theoretical restriction uiy /mz' < 0.26 



121. It turns out to be that this peculiar structure of the model imposes the 
25|, whereas lower bounds on my and mz' obtained from experimental data 



restrict this mass ratio to be 0.19 < niY /rnz' [25|. The behavior of the spin 1 amplitude as a function of the ray Imz' 
is shown in the third graph of Fig. [3] within the allowed interval 0.19 < my /mz' < 0.26. From this graph, it can be 
appreciated that |Mp ranges from 389 to 155 within this interval. At this level of amplitude, the vector contribution 
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is larger by more than two and three orders of magnitude than the fermionic and scalar contributions, respectively. 
As we will see below, this dominant effect of the vector particle is reinforced by the exotic charge contained of one of 
the bileptons. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



rrip /m^' /m^' 




150 ' ^ 

0.19 0.2 0.21 0.22 0.23 0.24 0.25 0.26 

FIG. 3: Sensitivity of the Z' — >■ 777 decay to the spin of the particle circulating in the loop. The combined effects of spin and 
mass ratio m/mz' are shown: first graphic, spin 1/2 effect, second graphic, spin effect, and third graphic, spin 1 effect. 



Having discussed the relative importance of the diverse spin contributions to the Z' — > 777 decay, we turn to 
discuss the corresponding branching ratio. To determine it, we need the main decay widths of the Z' predicted by the 
minimal 331 model, which are given in Ref. [l7j . As already commented, a good strategy for studying the sensitivity 
of this decay to physics lying at the mz> scale consists in assuming a degenerate spectrum for each type of new 
particles, i.e., we will assume that the three new exotic quarks are mass degenerate, the same will be assumed for the 
case of the four simply charged and three doubly charged Higgs bosons. The pairs of bileptons and also will 
be taken with the same mass. This assumption is quite reasonable, as all these new particles receive their mass at 
the firs stage of spontaneous symmetry breaking. The contribution to the Z' — > 777 decay of the exotic quarks and 
the SM fermions are shown in Figs. 2] and [SI respectively. From these figures, it can be appreciated that the exotic 
quarks contribution ranges from 10^^ to 4 x 10~^ for 0.1 < mg/m^' < 0.5, which is about three orders of magnitude 
larger than the combined contribution of the SM fermions. This is a surprising result, as it seems contradicts our 
previous analysis concerning the spin 1/2 contribution to this decay, in which we concluded that the contributions of 
the lightest SM fermions dominate. This apparent contradiction is due to the fact that the global factor involving 
products of coupling constants (see Eq. ([M)) ) is much larger in the case of exotic quarks than in the case of SM 
fermions, as it is shown in Table IIVI 

We now turn to discuss the contribution of the charged scalars to the Br(Z' — > 777). Its behavior as a function of 
tuh/tuz' is shown in Fig. [S] From this figure, it can be appreciated that Br{Z' — ^ 777) is larger for lowest values 
of niH/mz', which is consistent with our previous analysis for spin amplitudes shown in Fig. [2l The contribution 
to the branching ratio can be of order of 10~^, at the best. This contribution to the Z' — > 777 decay is one order of 
magnitude lower than the exotic quarks one and, as it can be appreciated from Table |Vl the exotic doubly charged 
scalars dominate. 
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FIG. 4: Exotic quark contribution to Br{Z' — ^ 777) as a function of niQ/mz' 
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FIG. 5: Contribution of the SM fermions to Br{Z' — >■ 777) as a function of m^/. The left (right)-handed graph corresponds 
to a scenario with mq — 500 GeV (mg = 1000 GeV). The lepton, quark, interference, and total contributions are separately 
shown. 



TABLE IV: Relative importance of global factor given by Eq. (|36p . 

I 0.096 

u,c 0.221 

d,s 0.023 

D,S 441.715 

b 0.007 

t 5.475 

T 3435.560 



As far as the bilepton contribution to the Z' — )■ 777 decay is concerned, the behavior of the branching ratio as a 
function of my /mz' is sho"wn in Fig. [71 From this figure, it can be seen that the bilepton contribution to Br{Z' — > 777) 
ranges, in the variation for my /mz' allo"wed by theoretical and experimental constraints, from approximately 8 x 10~^ 
to 3x 10~^ (see Tables IVll and IVlIl) . at the best. This contribution is larger than the exotic quark one by approximately 
a factor of 2. As it occurs in the case of exotic quarks, the exotic charge content of one of the bileptons play a decisive 
role in obtaining this result, as it can be appreciated from Table IVlIIl 

Finally, the total contribution to Br{Z' — >■ 777) is displayed in Fig. [5] From this figure, it can be appreciated that 
the branching ratio ranges from 7.87 x 10"'' to 5.78 x 10~^ for 0.5 TeV< m'^ <3 TeV, in a scenario with mg — 500 
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FIG. 6: Contribution of charged scalars to the Br{Z' — ^ 777) as a function of ranl'mz' for mg = 500 GeV (left-handed) and 
mg = 1000 GeV (right-handed). 



TABLE V: Relative importance of global factor given by Eq. (|49|l . 



H 




ht,ht 


1.025 


h+ 





ht 


0.028 


dr 


130.435 


d++ 


167.956 


4+ 


65.608 




0.19 0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.19 0.2 0.21 0.22 0.23 0.24 0.25 0.26 



my /m^' my/rn^' 

FIG. 7: Contribution of bilepton gauge bosons to Br{Z' — >■ 777) as a function of my /m'z for mq = 500 GeV (left-handed) 
and mq = 1000 GeV (right-handed). 



TABLE VI: Some values of Br{Z' 777) in the scenario mq = 500 GeV. 







rriY/mz' 




Br 




0.19 


0.26 


mz' 


= 1 TeV 


7.68x10'^ 


3.05x10^^ 


mz' 


=2 TeV 


2.73x10"^ 


1.09x10"^ 


rriz' 


=3 TeV 


2.61x10"^ 


1.04x10"^ 
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TABLE VII: Some values of Br{Z' 777) in the scenario mg = 1000 GeV. 







my /mzi 




Br 




0.19 


0.26 


mzi 


=1 TeV 


7.68x10"^ 


3.06x10"^ 


mzi 


=2 TeV 


7.67x10"'^ 


3.05x10^^ 


mz' 


=3 TeV 


2.97x10"^ 


1.18x10"^ 



TABLE VIII: Relative importance of global factor given by Eq. (|4T|) . 

V {gYf 

Y+ 0.226 
F++ 14.469 



GeV, my jm'^^ — 0.19, and mn — 250 GeV. As it is shown in Table HXl inside of this range of variation of mz' the 
branching ratio can reach a maximum of 1.07 x 10~^, and a minimum of 2.03 x 10~^ value, which occur for = 1.45 
TeV and 1 TeV, respectively. However, it can be appreciated from Fig. [5] that with the exception of a small interval 
centered in mz' = 1 TeV, the branching ratio for the Z' — > 777 decay is essentially of 10~^. It is worth comparing 
this result with the branching ratios obtained in reference [13] for the rare one- loop decays Z' — ^ ZZ and Z' Z7, 
which are of the order of 10~^ and 10~^°, respectively. This shows that the Z' — ^ ZZ and Z' — >■ 777 decays have 
branching ratios of the same order of magnitude in the minimal 331 model. This surprising result can be explained 
by noting that the Z' — > ZZ decay only receives contributions from the fermionic sector, whereas the Z' — >■ 777 one 
receives contributions from both the fermionic and bosonic sectors of the model. In addition, as already seen, the 
exotic charge content of the new quarks and gauge bosons play a central role in the three body decay. 

1.2x10'*' 

1.0x10"*' 

8.0x10'^ 

^ 6.0x10'^ 
m 

4.0x10'^ 
2.0x10'^ 
0.0x10° 

0.5 1 1.5 2 2.5 3 

m^. [TeV] 

FIG. 8: Total contribution to Br{Z' 777) for mg = 500 GeV, niY/mz' = 0.19 and ruH = 250 GeV. 



TABLE IX: Some values of the total Br(Z' 777) for mg = 500 GeV, mY/mz' = 0.19, and mn = 250 GeV (see Fig. 
The maximum value, 1.07 x 10~®, occurs at niz' ~ 1.45 TeV. 





rnz' [TcV] 




Br 


0.5 1 


3 


Total 


7.87x10"^ 2.03x10-'' 


5.78x10-'' 
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V. CONCLUSIONS 



Purely photonic decays of self-conjugate vector fields, to which the particle coincides with its antiparticle, as it is the 
case of the 1^" = Z, Z' gauge bosons considered in this work, are very constrained by gauge invariance, Bose statistics, 
and Lorentz invariance. Gauge invariance forbids any coupling of V'^ to photons at the tree level, so they only can 
arise as a quantum fluctuation of one-loop or higher orders. Since couplings of to two photons cannot exist due to 
the Landau- Yang theorem (26j , the interaction with three photons is the most important electromagnetic coupling of 
1^", which, in the context of a renormalizable theory, first arise at the one- loop level. Gauge invariance restricts this 
coupling to be characterized by dimension-six operators of the form {j j nx^yV^^^^F^" F\pF^f [l^l, where / represents 
a loop amplitude and m is the mass of the particle circulating in the loop. From the decoupling theorem [28|, one 
expects that if m 3> myo, the loop effect of the heavy particle decouples quickly. In contrast, one expect a relevant 
contribution if m < myo and it is more and more important if m is smaller than myo . In this work, we have studied 
the decays of Z and Z' into three photons within the context of the minimal 331 model, which predicts the existence of 
three new exotic quarks, two new gauge bosons, one simply charged and one doubly charged, and four simply charged 
and three doubly charged scalars. All the features of the ^"777 coupling commented above were reproduced. In the 
case of the SM Z gauge boson, it was found that the Z — )• 777 decay is insensitive to new physics effects, as the masses 
of the new particles are much larger than mz- Although large global factors arising from exotic charge contents of 
the new particles can substantially increase the loop amplitude, the well-known branching fraction of about 10~^, 
which is determined essentially by the lightest fermions and remains unchanged. Due to its insensitivity to heavy 
physics effects, this decay will likely be beyond the reach of the LHC or the future ILC. As far as the Z' — > 777 
decay is concerned, it was found that it can have a branching fraction as large as 10~^, which may be at the reach 
of future colliders. In particular, in the case of the LHC it has been found [1,] that the primary discovery mode for 
a Z' boson is a dilepton resonance via the Drell-Yan production process with a branching fraction of order 10^^ for 
a Z' mass in the range of the TeV scale and with an integrated luminosity of 100 to 300 ^h~^ . Our results for the 
branching fraction for Z' 777 thus makes it rather difficult to detect this decay mode at the LHC. In accordance 
with the previous general discussion, this decay is more important in the measure that the masses of the particles 
circulating in the loop turn out to be smaller than mz' ■ It was found that the exotic quarks and bilepton gauge bosons 
contribution is the dominant one, whereas the scalar contribution is smaller by about one order of magnitude. It was 
found that the exotic charge content of the particles circulating in the loop play a crucial role in this decay. In general 
terms, our study shows us that the decay of a self-conjugate vector boson V'^ into three photons is favored by three 
circumstances: 1) the presence of particles with exotic content of charge, 2) such particles have masses substantially 
smaller than wyo, and 3) there are several exotic particles. 
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Appendix A: Feynman rules 



In this Appendix, we present the Feynman rules used in the paper. 

In FigMg-iBiB = -ieQB, 9v°b^b = -i99v» I'^^^w , g^-^B^B = -^e^Ql, g^yoB^B = "ieQBggvof^cw, g-yS^S = 'i^Qs, 
9v°S'>S — *55yo/2cvi/, 577Sts = and gjv"sis — i^Qsggy'^/'^cw- On the other hand, the tensorial functions 

are given by 

rQiQ2a3(^l' ^2, fcs) = (^2 — ^3)qi5q2Q3 + ^ ~ ^1 ^ ^2 + ^3^ ffaiQa 

fcl - fc2 + "^^^-S fcg j g^^^^ ^ (Al) 

3*5 _ o„ „ (11 ^iNB'^B 

" aia2ti3a4 

T^^'^{k2,k:i) = {k2-h)o., (A3) 
where azw-w+ = 0i-yzw-w+ — oi-jB^b = ot-f-yB^B — ctv"Y^Y — a^yoyty — —1. 



010203 a4 — ^5010250304 I ^+ f K^o 103^02 04 + 3010450203 j I l^^j 
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Nr. 



BUk2 





FIG. 9: Bosonic Feynman rules needed for the calculation of the — >■ 777 decay. In this figures, AT = 7, V°, with V° 

Z, Z'; B = W^ , , y*"^; S = G^, C^, C^, Gy, Cy, Cy, Gy^, Cy"*", Cy'^ S = H = h^ , ftj, /ig", ft^ , rf]^'^, rf^'^, 



The diverse factors appearing in Sec. IIIII concerning the coupHngs of the Z' gauge boson to pairs of fermions, 
charges scalars, and bileptons are given in tables [Xj |Xll and Kill respectively. 



TABLE X: Couplings of Z and Z' gauge bosons to pairs of SM and exotic fermions. Only the vectorial component is shown. 



/ 


Qf 


9vz 




e,n,T 
u, c 


-1 

2 


2 

3-83^ 


2c^ 


3 


6 




d, s 
D,S 
b 
t 

T 


1 


3-4sf,, 


1 


3 

4 
3 

1 


6 

3 

3-4sf,, 




3 

2 


6 


l+4s^ 


3 
5 


6 




3 


3 





TABLE XI: Couplings of Z and Z' gauge bosons to pairs of charged scalars. The value S ~ is assumed for the Z' ~ Z mixing 
angle. Also the values jAij, j/jj ~ 1 are assumed. 



s 


Qs 


gl 


9%. 


G%f 


1 






Gy 


1 


— 1 — 2SiY 


^{l-2tl) 




2 




£^(l-2ti) 


ht,h+ 


1 


— 2s^ 


2cn/ce 


h+ 


1 


C21V 


^-^[l-2{l + tl)HU 


ht 


1 


C21V 






2 


-4s^ 


2c^c,^2j2x| + (l tlYa^] 


4+ 


2 




2^iV|[2X| + (l-t^)a2] 


4+ 


2 


2C2W 


2c\vcf) 
^/3 
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TABLE XII: Couplings of Z and Z' to pairs of charged gauge bosons. 



B 


Qb 


fff 


5f 




W+ 


1 









Y+ 


1 


-1 - 2slv 






Y++ 


2 


1 - 44/ 


%/3\/l- 





Appendix B: Form factors 

The form factors associated with the spin 1/2 contributions are given in Ref . [l^ . Here we present the form factors 
induced by the bosonic particles of the model. The spin 1 form factors are given by 

^ g(f(l)(pi3-2p23)a ^ ■Bo''(2)(pi2+P23)« g(f(3)(2pi2-pi3)p23a 
12p?3P23 12j5?2P23 12p22p2g 

■^^(4)(Pl2 +Pl3 +P23)[Pl2(Pl3 - 2p23) +P13P23]Q 
12p?2Pl3?'23 

, Co^(l){[abi3 - 2pi3) + 12pi3Pl3]P?2 + 3m|p2gp2^a| C^f (2)pi3 [a(p?2 + P23) - 12pi3Pi2] 



12pi2P?3P23 12P?2P23 
^ C^(3)p23[2(6pi3 - P23a)Pi2 + 3pj3(m|Q! - 4pi3)pi2 +Pi3P23a] 

12P?2P?3 

g(f (4)bi3 +P23)(apj3 + 12pi3Pi3 - 2p33a) 
12Pl3Pi3 

C^(5)(pi2 +P23)[«(Pi3P?2 +Pi3Pi3) + 3pi3(m|a - 4^13)^12] 

12pi2Pl3Pi3 

C(f (6)(pi2 +Pl3)[2(6pi3 -P23a)pj2 - 12pf3Pi2 +Pi3P23a] 
12Pl2Pl3 

-D^(l)[2pi2P23(12p23 +Pi3a)™| +pi3(2m|Q - 3pi3/3)to| + 2pf2Pi3«] 

12pi2Pl3Pl3 

-Po^(2){2pi3aTO| +Pi3[4pi2(6pi3 -P23) + 3pi3(p23'a - 8pi3)jw| + 4^^2^23(6^13 "P23'a)} 

12pi2P?3 

-Do^(3)[2pi2am|; +pi2(24p^2 - 24pi3pi2 + hpi^p2za)m% + 2^^3^23(^230 - 12pi2)] 

12p?2Pl3 

■ 

^ B^(l)a g^(2)[pi2(2pi3+p23)+P23(3pi3+P23)]« (4) (pi2 + Pl3 + P23) (2pi2 + 3p23)« 

2p33 4p33(pi2+P23)2 4p33(pi2 + P23)2 

j '^^(1)(2P13 +P23)(P23W| +Pl2Pl3)« 
4P13P!3 

, C'^(2)bl3Pi3™B +Pl2P23(2pi3 +P23)m| +p?2Pl3(2pi3 +P23)]a , C^(3)m|a 



C^(4)(pi3 +P23)[2p23»7l| +Pl2(2pi3 +P23)]a 



4pj2P23 4pf2P23 

B 



4pi2P23 

Co^(5)[(2pi3 +P23)pf2 + 2p23(™| + 2pi3 +P23)pj2 +Pi3(3?Ttlj + 2pi3 +P23)Pi2 + m%p^^]a 

4P?2P23(P12 +P23) 

g^(6)m|(pi2 +P13)^« _ £'^(l)[2pi3"^B +P12P23(8P13 + 3p23)m| + 2pf2Pi3(2pi3 + p23)]a 
4p?2Pl3Pi3 4pi2p|3 

D§ {2)m%{2pi3m% +pi2P23)a I?^(3)m|[2pi2m| + (pi2 + 2pi3)p23]a a 



4pi2Pl3Pi3 4pf2Pi3 4p^3(pi2+P23 



(B2) 



S13 



+ 



Bp {l)[8pi2Pl3 + 3(4pi2 + Pl3)P23Pl3 - 6^12^13^13 - (4pi2 + ^Pl3)P23h 

24p?3?5i3(Pl3 +P23)^ 
B^{2){8pi2-P23)a go^(3)(4pi2 +Pl3)ct 

24pi2pl3 24pi2pf3 
-So^(4)(Pl2 +P13 +P23)[2pi2(4pi3 + 6p23Pl3 " 3^23^13 - 2p23) - Pl3P23(Pl3 + P23)^]Q: 

24pi2P?3P23(Pl3 +P23)^ 

C^(l)pi2{4pi2(2pf3 -pi3)a + 3pi3P23[a(2pi3m^ +^^3) -pia/?]} 

n3 „3 



24P?3P23 



Co^(2)[8pi3P?2 + 3(2m| +pi3)p23P?2 -P13P23]" C-q^ (3)p23 [a(4p?2 + P13) + 3P13^P12] 



24P?2P23 24p22p3g 
+ 0, 3 rV^' , 7{3?3l3P23[-2pi3(P?3 + 2^23^13 " ^13)""^! " (Pl3 + P23)^(P?3" " f'23/3)] 



go^(4) 
24Pi3P23(Pl3 +P23) 



+ 



-4pi2(]5l3 +P23)^(2p?3 -P^3)a} 
Co^(5)(Pl2 +P23)[Pl3Pl2 + 3(2m| +pi3)p23Pl2 - Pl3P23]« 
24p?2Pl3P23 

Co''(6)(pi2 +Pi3)[a(4p?2 +P13) + 3pi3/3p?2] 



24p?2P?a 
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■D^(l){a[8pj2Pl3 +Pi2(14m| + 3pi3)p23Pl3] + pi3[2aTO|; + 3(pi3Ct - 8p23)ml + 12pi3P23]} 

12P13P23 

+ 77^^ [-2p?3aTO| +Pl3P23(4pi2a + 3pi3/3)m| +P23(-12p?3 + 3pi2P23/3pi3 + 4pi2P23a)] 

l2pj3P23 

-Do'(3)[-2Pi2Q;m| +Pl2P23(24pi2 + J3l3a)m|5 +Pl3P23(Pl3P23Q: - 12pf2)] 



+ 



12pi2Pl3P23 

(2pi3 -P23)Q! 



12pi3P23(Pl3 +P23) 

In the above expressions, B stands for W~^, or Y'++. In addition, a = as — 3, /3 = as + 5, aw=t'y/, and 
On the other hand, the form factors associated with spin particles are given by 

^ ■B^(l)(j7i3-2p23) -Bo^(2)(pi2+P23) ^ -Bo^(3)(2pi2-pi3)P23 
12p?3P23 12p?2P23 12pi2P? 



+ 



3 

-Bp (4)(pi2 +P13 +P23)[P12(P13 - 2p23) + P13P23] Cp (l)(Pl2Pl3 + 3m|p|3P^3 - 2pf2P23) 



12pi2H3P23 12pi2Pi3P23 
C^(2)Pl3(Pl2 +P23) ^ C'^(3)P23(2p23Pl2 " 3m|pf3Pi2 - P13P23) 

12Pl2Pi3 12Pl2Pl3 
I Cp (4)(Pl3 +P23)(pj3 - 2P23) ^ gp (5)(Pl2 +P23)(Pl3P?2 + 3w|pi3Pi2 +P13P23) 
12Pl3Pi3 12p?2Pl3Pi3 
Co^(6)(pi2 +Pl3)(2p?2 -F?3)P23 _^ -D^ (1) [2^23^4 +pi3(2pi2 - 3^23)^23^4 + 2p?2P?3] 



+ 



12p?2Pl3 12pi2Pl3P23 

£)f (2)[2pf3m| +pi3(3pi3 - 4pi2)p23m| - 4pf2P23] 

12pi2P?3 
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^ -Dp (3)(2pi2m| +Pi3P23)(P12to| + 2pi3P23) 1^ 



12P12P13 12pi2Pl3 
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f ^^(1) g^(2)[pi2(2pi3+P23)+P23(3pi3+P23)] g^(4)(pi2 + Pl3 + P23)(2pi2 + 3^23) 

'2ph ^pUPl2+P23r 4p33(pi2+P23)2 

Cq {l){2pi3 + P23){p23ml + P12P13) 

'ipwpis 

C^(2){m|[pi3P23 +Pi2P23(2pi3 +P23)] +Pi2Pi3(2pi3 +P23)} _ C^(3)m| 

4p?2P23 4p?2P23 
Co'(4)(pi3 +P23)[2P23W| +Pi2(2pi3 +P23)] 
4j3l2P|3 

C^(5)[(2pi3 +p23)Pi2 + 2p23(w| + 2pi3 + P23) P 12 + phi^m^ + 2pi3 +P23)Pi2 + "^1^23] 

4P?2^'23(^'12 +P23) 



I go'(6)m|(pi2+Pl3)^ ^ £'f(l)[2pi3m|+pi2P23(8pi3 + 3p23)m| + 2p^2Pl3(2pi3+P23)] 

4p?2Pl3Pi3 4pi2]5t3 

^ I^g(2)(2pi3m|+pi2P23w|) ^ j^^X3)m|[2pi2m| + (pi2 + 2^13)^23] ^ 1 

4pi2Pl3Pi3 4pf2Pi3 4p|3(pi2 +J)23) ' 



/S13 = 



^0 (l)bl2(6p23Pl3 - 8P13) - 3(4pi2 +P13)P23P13 + (4pi2 + 3pi3)pi3] 

24p?3P23(Pl3 +P23)^ 
^ Eo^(2)(p23-8pi2) ^ go^(3)(4pi2+Pl3) 
24pi2Pi3 24pi2P?3 
-Bo(4)(pi2 +P13 +P23)[2pi2(4pf3 + 6p23Pi3 " 3^13^13 - 2^13) - pi3P23(pi3 + P23)^] 

24pi2pf3Pi3(pi3 +P23)^ 
<^^(1)P12[-8P12P13 - 3(2m| + Pl3)P23Pl3 + (4pi2 + 3pi3)P23] 

24p?3l'i3 

^ C^{2)[-8pi3pl2 - 3(2m| +pi3)p23Pi2 +P13P23] I go^(3)(4pf2 + 3pi3Pi2 + ^^3)^23 
24pf2Pi3 24p?2Pl3 

+ 0, 3 3 V , r{4pi2(2p?3 -P^3)(pi3 +P23)^ +3pi3P23[(Pl3 -P23)(P13 +P23)^ 



2 /2 ,0 2 Ml , <^0(5)(Pl2+P23)bl3(8Pi2~P23) + 3(2m|+pi3)p23Pl2] 



+2m,Pi3(Pl3 + 2P23P13 - P23)]} + ■ 24p^2Pl3P^3 
Co^(6)(Pl2 +Pl3)^(4p?2 -P13P12 +P13) 



+ 



24p?2P?3 

^0 (l)[8Pi2Pi3 +Pi2(14m| + 3pi3)p23Pi3 + m|(2m| + 3pi3)pi3] 

12pi3P23 

£>g(2)[-2pf3m| +Pi3(4pi2 + 3pi3)p23m| +Pl2(4pi2 + 3pi3)pi3] 

n3 



12Pi3P23 

I -Dp (3)(m|pi2 - Pi3P23)(2pi2m| +P13P23) _ P23 - 2pi3 .gg. 

12pf2Pl3P23 12pi3P23(Pl3 +P23) ' 

where S—h^, /ij, /ij, /14 , rf^"^, dj"'', £^3 ^, G^, C^, C^, Gy, Cy, Cy, Gy^ , Gy^ , Cy"*". In the above expressions, 
we have introduced the following defintions for the Passarino-Veltman scalar functions: 

B^{1) = Bo{2pt2,m\,m\), B^{2) = Bo{2pr3,m\,m\), B^{3) = Bo{2p23,m\,m\), 
B^{A) = Bo{mlo,m\,m\), G^ {1) = Go{Q,Q,2pi2,m\,m\,m\), G^ {2) = Co(0, 0, 2pi3, m^, m^, m^), 
(3) = Co (0, , 2p23 , , m| , m| ) , (4) = Co (0 , 2pi2 , m2,o , , m| , ) , (5) = Co (0 , 2pi3 , m2,o , m| , , ) , 
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C^{&) = Co{0,2p23,my,mx,mx,mj^), (!) = i?o(0, 0, 0, m^o, 2pi2, 2pi3, m^, m^, m^, to^); 

^cf (2) = Do{0,0,0,myo,2pi2,2p23,mx,mx,mx,mx), (3) = Do(0,0,0,myo,2pi3,2p23,mx,mx,mx,mx), 
where X denotes the virtual particle circulating in the loop, and p^j = pi ■ pj with i,j=l, 2, 3. 
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